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Abstract 

We investigate the spin polarization of r ± leptons produced in v T and 
v T nucleon scattering via charged currents. Quasi-elastic scattering, A res- 
onance production and deep inelastic scattering processes are studied. The 
polarization information is essential for measuring the appearance rate 
in long baseline neutrino oscillation experiments, because the decay particle 
distributions depend crucially on the r 1 * 1 spin. In this article, we calculate 
the spin density matrix of each process and estimate the spin polarization 
vector in medium and high neutrino energy interactions. We find that the 
produced r^s have high degree of polarization, and their spin direction 
depends non-trivially on the energy and the scattering angle of in the 
laboratory frame. 
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1 Introduction 



Recent studies from neutrino oscillation experiments are revealing the amazing 
nature of the neutrino sector, with their non-zero masses and large flavor mixings. Es- 
pecially, reports from Super-Kamiokande (SK) collaboration [I] strongly suggest that 
nearly maximal oscillation from into v T is occurring in the atmospheric neutrino 
flux. To demonstrate this oscillation, it is important to detect v T appearance in os- 
cillation experiments. Several long-baseline neutrino oscillation experiments, such as 
ICARUS [2], MINOS 0, OPERA @] are proposed, and they are expected to detect the 
r appearance by charged current (CC) reactions off a nucleon target 



with N = p,n. Because r production by a nucleon target has a threshold for neutrino 
energy at E v rs 3.5GeV, these experiments should provide high energy neutrino flux. 
It has also been pointed out by Hall and Murayamaj^j, that SK may be able to detect 
the t appearance events with more than several years of running. 

The produced r decays into several particles, always including a neutrino (v T ). 
Therefore the r appearance signal should be obtained from decay particle distributions. 
Because the r decay distributions depend significantly on it's spin polarization the 
polarization information is essential for us to identify the r production signal. 

r polarization should also be studied in order to estimate background events for 
the Up — > v e appearance reaction, which will be searched for in neutrino oscillation 
experiments, such as those using high intensity neutrino beams from J-PARC[7J. Be- 
cause the oscillation amplitude of — > v T is larger than that of — > t'ejHj, and 
because the branching ratio of r T — > e T + X is relatively large, the e production via 
the — > v T — > r — > e chain can be significant [H] . Since the e energy and angular 
distribution depends crucially on the r polarization, it's information is necessary to 
estimate the background. 

So far, several authors have calculated the r production cross section for nucleon 
targets |51 ITUl ITTj . but to our knowledge, no estimation of the polarization of produced 
r's is available. In this paper, we study the spin polarization of r produced by v T scat- 
tering off a nucleon target. We consider the quasi-elastic scattering (QE), A resonance 
production (RES), and deep inelastic scattering (DIS) processes, which are known to 
give dominant contributions in the medium and high neutrino energy region [TO]. The 
spin polarization vector is obtained from the spin density matrix which is calculated 
for each process. 

The article is organized as follows. We give the general kinematics of r production 
in neutrino-nucleon interaction and the relation between the spin density matrix and 
the spin polarization vector in section |21 Then we present the details of the spin density 




(1) 
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matrix calculation of QE, RES, and DIS processes, in sections El IH andEl respectively. 
In section El the differential cross section and the spin polarization vector of produced 
are estimated for medium and high neutrino energies. Section |7| gives discussions 
and our conclusions. 



2 Kinematics and Formalism 



In this section, we show the physical regions of kinematical variables and give the 
relation of the r spin polarization vector and the spin density matrix of the charged 
current (CC) r production process. Firstly, we define the four-momenta of incoming 
neutrino (k), target nucleon (p) and produced r lepton (k') in the laboratory frame 



F = {E v ,0,0,E v 
= (M, 0,0,0), 



(2) 
(3) 
(4) 



k'^ = (E T ,p T sm9,0,p T cos8). 

Here, E v and E T are the incoming neutrino and outgoing r energies, respectively, in 
the laboratory frame, M is the nucleon mass, and p T = yj ' E 2 — m 2 with the r lepton 
mass m T = 1.78GeV. We also define some Lorentz invariant variables 



Q 2 = -g 2 , q» = k»- fc'", 
W 2 = (p + qf. 



(5) 
(6) 



Q 2 is the magnitude of the momentum transfer and W is the hadronic invariant mass. 
The physical regions of these variables are given by 

M < W < yfl -m T , 



(7) 



and 



where s = (k + p) 2 and 



Q 2 .(W)<Q 2 <Q 2 AW), 



Qi(W) 



- M 2 1 

— (1±P) - - 
2 v H! 2 



W 2 + m 2_M_(w^ m ^ 



(9) 



with [} = A5 (1, m 2 /s, W 2 /s) and X(a, b, c) = a 2 + b 2 + c 2 - 2(ab + bc + ca). 

The scaling variables are defined as usual: 

Q 2 Q 2 



X 



y 



2p-q W 2 + Q 2 -M 2 ' 
W 2 + Q 2 - M 2 



P-Q 
p-k 



s-M 2 
2 



-I- 



(10) 



Here, x is the Bjorken variable and y is the inelasticy. The physical regions for x and 
y are obtained by Albright and Jarlskog |lH I12j: 



mz 



2M{E V - m 7 



< x < 1 



(12) 



and 



where 



A- B <y < A + B, 



A=l(i- T - £± 
2 I 2ME v i 2£^ 



B= l - 
2 



+ 



2 2 



2ME v x J El 



1 + 



2E, 



(13) 

(14) 
(15) 



The above regions agree with those determined by Eq.((Zj) and Eq.(jEJ). 

We label the relevant subprocesses by using the hadronic invariant mass W and the 
momentum transfer Q 2 . We label QE (quasi-elastic scattering) when W = M, RES 
(resonance production) when M + m n < W < W CVLt , and IS (inelastic scattering) when 
W cut < W < ^fs — m T . W cvX is an artificial boundary between RES and IS processes, to 
avoid double counting. The W cv x value is taken in the region 1.4GeV~1.6GeV. Within 
the IS region, the region where Q 2 > lGeV 2 may be labeled as DIS, where the use of 
the parton model can be justified. 



FigH shows the kinematical regions of each QE, RES and IS process on the x-y 
plane (left) and the p T cos 8-p T sin 9 plane (right) at E v = 10 GeV. The QE region is 
shown by open circles, the RES region by open triangles, and the DIS region is shown 
by the cross symbols. The region shown by the star symbol (*) gives the IS process 
at low Q 2 (Q 2 < lGeV 2 ). In this region the parton model is not reliable and we must 
use the experimental data to reduce errors. In this report, however, we use the parton 
model throughout the IS region. Studies on uncertainties in this region will be reported 
elsewhere. 
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Figure 1: Physical region at E u = 10 GeV in the x-y plane (left), and in the p T cosO- 
p T sin 9 plane (right). Open circles denote QE (quasi-elastic scattering), open triangles 
denote RES (A resonance production), and the cross symbols give the DIS (deep 
inelastic scattering) region with Q 2 > lGeV 2 . The region marked by the star symbol 
(*) gives inelastic scattering (IS), at W > 1.4GeV and Q 2 < lGeV 2 . 



Produced r will be partially polarized. We define the spin polarization vector 
parameterized as 

, P 



\Sx 1 Syi S z i 



(sin 6p cos (fp, sin dp simp P , cos6>p) 



(16) 



in the r rest frame in which the z-axis is taken along it's momentum direction in the 
laboratory frame. In Eq. (jl6|) . 9p and (fp are the polar and azimuthal angle of the 
spin vector in the r rest frame, respectively, and P denotes the degree of polarization. 
P = 1 gives the fully polarized r, and P = gives unpolarized r. The azimuthal angle 
is measured from the scattering plane where tpp = n/2 is along the p u x p T direction in 
the laboratory frame. The degree of polarization (P) and the spin directions (9p,<pp) 
are functions of E T and cos#. This spin polarization vector is related with the spin 
density matrix R\y, by the following relation: 

dRxw If 1 + PcosOp P sin Op e ilpp \ dcr snm 
dE T dcosO 2\ Pwn9 P e-*e p 1 - P cos0 P ) ' dE T dcosO' 1 ' 

The density matrix is calculated as R\y oc J2M\M y , where M\ is the helicity ampli- 
tude with the t helicity A/2 defined in the laboratory frame, and da sum = dR ++ + dR 

is the usual spin summed cross section. The summation symbol implies the summation 
over final states, and the spins of the target and final-state particles. 



The spin density matrix of r production is obtained by using the leptonic and 
hadronic tensor as 

dRxy G 2 f k 2 p T 



dE T d cos 



T^ v W 
4tt ME V ' 



where Gf is Fermi constant and k = M^/i^Q 2 + M^) is the propagator factor with 
the P^-boson mass My/ = 80.4GeV. For r~ production, the leptonic tensor V^ x , is 
expressed as 

LTy=fxJ7, (19) 
where the leptonic weak current j£ is 



2E U (E T -p T ) (sinf,- cosf, i cosf, sinf) (A = +) 
2E V {E T + p T ) (cosf, sin|, —i sin|, cosf) (A = -) 

in the laboratory frame. For r + production, we must replace the leptonic tensor L^, 
into defined as 

LZ> (21) 

where j£ is 

Tx=v v {k) 1 »' ] -^v T {k\\) 



2E u (E T +p T ) (cosf , sinf , i sinf, cosf ) (A = +) 

2E U (E T - p T ) (— sin|, cosf, % cosf, —sinf) (A = -), 

which is related with j£ by j£ = Aj*^, in the phase convention of Ref.[13j. In the 
following sections, we will abbreviate the overline of the leptonic tensor and currents 
for r + production process. 

The hadronic tensor is expressed in general as 

WV(P, ?) = -g^ip-q, Q 2 ) + ^ W 2 (p-q, Q 2 ) - ie^ap^ W 3 (p-q, Q 2 ) 

+^ W<(p.q, Q 2 ) + m \ + M q f V W 5 (p.q, Q% (23) 

where the totally anti-symmetric tensor e^ va p is defined as e i23 = 1, and the structure 
functions Wj=i v .. i5 (p • q, Q 2 ) can be estimated for each subprocess. Since q ■ j\ is pro- 
portional to m T , the structure functions W4 and W 5 appear only in the heavy lepton 
production case [12]. 



Inserting these equations into Eq.(fTH|) and Eq.(JT7J), we find 

dOsum G 2 f k 2 p T \ ( ^ JT t m 2 

dE T dcos9 2tt M 



± 

G 2 f k 2 p T 



{ ( 2Tyl + Wi ) (Er " Pr cos e) + W2 (Er + Pr cos 6) 

\Al 771 ^ ^ 



2. M- < M » 



and the spin polarization vector takes 

„ = T Vhfl ( 2Wl _ w , ± !k ^ - ^ + | /F, (25a) 
s y = 0, (25b) 
= T ~ | (2 W x - ^ Wi) ( Pr - E T cos 6) + W 2 (p T + E T cos 9) 

1~ ({E v + E T )p T - (E V E T + pi) cos 9) - ^ W 5 cos / F, (25c) 
for t t productions. The degree of polarization is given by 

p = V s ' + a l + «2- ( 26 ) 

The above results Eq.tj24 j) -tj26j ) agree with Ref.|12j. From the above equations we find 
(i) (fp takes either or tt for any r momentum, which means that the polarization 
vector lies always in the scattering plane, and (ii) if m T = 0, then s could take only 
(0,0, =f±), which means fully left-handed r or right-handed r + . 



3 Quasi-Elastic Scattering 



In this section, we give the spin density matrix calculation for the QE scattering 
processes 



v T + n —y t + p 
v T +p -> 



r + + n. 



(27) 
(28) 



Following Llewellyn Smith [14], the hadronic tensor is written by using the weak tran- 
sition current J,^ as follows: 



hi/ 



4 z — ' f 

spins 



(29) 



where # c is the Cabibbo angle. The weak transition currents J / ^ + ' ) and J( ^ for the iv 
and v T scattering, respectively, are defined as 

■J { , +) = W)\J { ^n{p))=u p (p')T^p\p)u n {p), (30) 

JH = (n(?/)| jjf>|p(p)) = SnCp'Jr^^tipCp) = (p(p)|j£ +) WpT, (31) 

where Y ^ is written in terms of the six weak form factors of the nucleon, F^ 2 3, Fa, F% 
and F p , as 



+ 



7^A(g 2 ) + 



(p + p'] 

M 



75- 



(32) 
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For the v T scattering, the vertex is obtained by r M (p',p) = •y Q Tj 1 (p,p')j . We can 
drop two form factors, FY and F^, because of time reversal invariance and isospin 
symmetry (or equivalently no second-class currents). Moreover, the vector form factor 
FY and FY are related to the electromagnetic form factors of nucleons under the 
conserved vector current (CVC) hypothesis: 

a 2 

G v F (q 2 ) ~ ^—GXJq 2 ) r v ir?\ r v (r?\ 

*{V)= m l iFZif) = G *^-%*to\ (33) 

1 - 1 - ^ 



AM 2 AM 2 



where 



GW) = ~ ( j— 2, GlM Z ) = - ( '-±y> (34) 

I 1 " M) \ m) 

with a vector mass My = 0.84 GeV and £ = /x p — /i n = 3.706. /x p and /i n are the 

anomalous magnetic moments of proton and neutron, respectively. For the axial vector 
form factor Fa, we adopt the following parametrization: 

F A tf) = . (35) 
1- q ^ 



with an axial- vector mass Ma = 1.0 GeV and Fa{0) = — 1.23 [I^j. For the pseudo-scalar 
form factor F p , we adopt the parametrization of Ref.f 



F p (g 2 ) = 2M 2 ^1 (36) 

with the pion mass m n = 0.14 GeV. The normalization of F p (0) is fixed by the partially 
conserved axial vector current (PCAC) hypothesis. It should be stressed here that, the 
form factor F p (q 2 ) has not been measured experimentally because its contribution is 
proportional to the lepton mass*. The production cross section and the polarization 
of r are sensitive to F p (q 2 ) because of the large r mass and the spin-flip nature of the 
form factor. 



In FigEJ we show the total cross sections of the QE process versus the incoming 
neutrino energy. We plot not only the r-neutrino interaction process, but also the /Li- 
neutrino interaction process for comparison. Solid curves are the v p and v T scattering 

*After the paper[22] was published, we learned that there are some experiments which measure the 
pseudoscalar form factors in muon capture 23 and in pion electroproduction(24| . Those experiments 
found results consistent with the PCAC relation at very low Q 2 , but they are not sensitive to the 
large Q 2 region (w 0(lGeV 2 )) which is relevant for t production. A lattice study by Liu et a/.|2"5] 
seems to agree with our parametrization of Eq. . 
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cross sections and dashed curves are the and v T scattering cross sections. Our results 
agree well with those of Hall and MurayamajS]. 



1.2 




Figure 2: The neutrino energy dependence of the total cross sections of the QE (quasi- 
elastic) processes, vin — > l~p (solid lines) and \>ip — > l + n (dashed lines). The thick lines 
are for I = r, while the thin lines are for I = fi. 



4 Resonance Production 

In this section, we present the spin density matrix calculation for the A production 
processes 

u T + n(p) -> T~ + A+(A++), (37) 
u T +p(n) ^r+ + A°(A-). (38) 

We neglect iV* and the other higher resonance states, which are known to give small 
contributions [II H H5 | UK]. For the A resonance production, we calculate the hadronic 
tensor by using the nucleon-A weak transition current J M as follows: 

cos 2 g c ^ I WT(W) 

v - 4 ^s " "-* (w*-Miy + w^{wy { ' 

Here, Ma is the resonance mass, Ma = 1.232GeV, and T(H / ) is its running width 
estimated by assuming the dominance of S-wave A — > + tt decay: 

with T(M A ) = 0.12 GeV and X(a, b, c) = a 2 + b 2 + c 2 - 2(ab + be + ca). 
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The current J M for the process v T + n — ► r + A + is defined by 



J, = (A + ( P ')\JMp)} = ^ip') T» a u n (p), 



(41) 



where ip a is the spin-3/2 particle wave function and the vertex is expressed in 
terms of the eight weak form factors CYj£ t s 6 [T H IT7 | ITH | fTHj as 



M 



iff mi ~ ln<la) + {QtmP'-q - P^q 

c v c v 



+ 



M 



7s 



- i^q a ) + [p m p'-q - p'^qc 



M 2 

+C5 9/ia + j-p q^q a - 



(42) 



By using the isospin invariance and the Wigner-Eckart theorem, we obtain the other 
nucleon-A weak transition currents as 



(A ++ |J» = V3(A+| j» = >/3(A°|J M b> = (A-| J». 

From the CVC hypothesis, C\ = and the other vector form factors C^ 3 4 5 
to the electromagnetic form factors. We adopt the following parametrizations: 



(43) 
are related 



c 3 v (o) 



2 • 



M 
Ma 



<W), <W) = o, 



(44) 



with 63(0) = 2.05 and a vector mass My = 0.54 GeV 2 . For the axial vector form 
factors Ql 345 , we use the modified dipole form factors [T7| ITH] 



'{=3,4,5 



= Ci(0) 



Mi, 



(45) 



with C 3 (0) = 0, C 4 (0) = -0.3, C 5 (0) = 1.2, a 4 = a 5 = -1.21, 64 
Ma = 1.0 GeV. And for Cg , we use the following relation [2*0]: 



65 = 2.0GeV 2 and 



OV) = Ct{q 2 ) 



M 2 



9 ; 

— q z 



(46) 



which agrees with the off-diagonal Goldberger-Treiman relation in the limit of m 2 — > 
and g 2 — > 0|TZ|. The pseudo-scalar form factor C§(q 2 ) has not been measured because 
its contribution vanishes for massless leptons. As in the case of the F p (q 2 ) form factor 
of the QE process, C^(g 2 ) has significant effects on the r production cross section and 
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the r polarization. 



In Eq. (j39J) . summation over the hadronic spins is done by using a spin projection 
operator of the spin-3/2 particle wave function which is given by 



spin 



2 p'f+p'y I p'^lv-p'vlv 

3 Ml 3 



The hadronic tensor is now calculated by 



W RES = C ° s2 ° c Tr 
V" 4 



u/3 



WT(W) 



Tr (w 2 - Ml) 2 + w 2 r 2 (wy 



(47) 



(4* 



By integrating over E T and cos 6 within the kinematical region of M + < W < 
1.4GeV, we estimate the total cross section of the A production (RES) processes. In 
Fig |2J we show the total cross section versus the incoming neutrino energy. We also 
plot the total A production cross sections for and scattering processes, in order to 
examine the lepton mass dependence. The ^ production cross sections grow sharply 
at low E v , while the r production cross section grow mildly from around E v = 4GeV. 
The cross sections of A ++ and A~ production processes are larger than those of A + and 
A productions. This feature is expected from the Clebsh-Gordan coefficients of the 
transition currents, in Eq. (j43|) . Our results agree approximately with those of Paschos 
and Yu[TU], which include the contributions from N*(Sn, P\\) resonance productions. 




Figure 3: Total cross sections of the A production (RES) processes, plotted against 
the incoming (anti) neutrino energy. The solid, long-dashed, dashed, and dot-dashed 
lines show A ++ , A + , A , and A~ production cross sections, respectively. The thick 
lines are for v T and v T scatterings and the thin lines are for and scatterings. 
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5 Deep Inelastic Scattering 



In this section, we present the spin density matrix calculation for the deep inelastic 
scattering (DIS) processes 

v T + iV -> T~ + X, (49) 
u T + N -> r+ + X. (50) 



In the DIS region, the hadronic tensor is estimated by using the quark-parton model; 

q,q 



W™(p, q) = £ / f f q ,,(t Q 2 ) KjtPfa, ?)■ (51) 



Here, p^ = £p M is the four-momentum of the scattering quark, £ is its momentum 
fraction, and f q and f q are the parton distribution function(PDF)'s inside a nucleon. 
By taking the spin average of initial quark and by summing over the final quark spins, 
we find the quark tensor 



K^(p q ,q) = 6(2p q -q-Q 2 -m 2 ) 



x2 



-9tx V {Pq-q) + tPq^Pqu =F ie^apPqqP + {Pq^v + q^Pqu) ■ (52) 



The upper sign should be taken for quarks and the lower for antiquarks. We retain 
the final quark mass, m q ', for the charm quark as m c = 1.5GeV, but otherwise we set 
m q i = 0. We neglect charm and heavier-quark distributions in the nucleon, as well as 
bottom and top production cross sections. 

By neglecting the nucleon mass and the initial quark masses consistently, we find 
the following relations: 

M 2 

W 1 (p-q,Q 2 )=F 1 (x,Q 2 ), W t=2 _ 5 (p-q,Q 2 ) = — - F i=2 _ 5 (x, Q 2 ). (53) 

Here, 

Fi = £/ M (£,Q 2 ), (54a) 

i,g 

F 2 = 2^Uq,q(^Q 2 ), (54b) 

f 3 = 2E Q 2 ) - 2 £ Q 2 )' (54c) 
<? g 

F A = , (54d) 
F 5 = 2£/ M -(£,Q 2 ), (54e) 

where the momentum fraction is £ = x for massless final quarks (m q i = 0), and £ = x/X 
with A = Q 2 /(Q 2 + m 2 ,) for q' = c. In the m c — > limit, the Callan-Gross relation 
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i<2 = 2xF\ and the Albright- Jarlskog relations = 0, 2xF$ = F2 hold. 

However, the differential cross section (Eq. (|24)l ) does not satisfy the positivity con- 
dition near the threshold with this naive replacement. We find that the following 
modification of the W\ structure function suffices to ensure the positivity constraints^: 

W x = ( 1 + ^) F lf (55) 

and find that the positivity is maintained when the charm quark mass is introduced 
by using the rescaling variable £ = x/X. 

There is further uncertainty in our parton model predictions for the inelastic scatter- 
ing processes where the hadronic final state is heavy, W > 1.4GeV, but the momentum 
transfer is small, Q 2 < lGeV 2 . This is the region of the phase space depicted by the 
star symbol (*) in the x-y plane and the p T cos 6-p T sin 9 plane of Figd In order to 
estimate the cross section and the spin polarization vector in this region, we use naive 
extrapolation of the parton model calculation, by using the parton distribution at the 
minimum Q 2 (Qjj = 1.25GeV 2 for the parametrization of A.D.Martin et a/. [21]) even 
when Q 2 < Ql- 



1 



/ v, + p -> / + X 

v ( + n -> r + X 




E, (GeV) 



Figure 4: Total cross sections of the DIS processes divided by the neutrino energy are 
plotted against the neutrino energy. Solid, long-dashed, dashed, and dot-dashed lines 
show Uip — > l~X, vin — > l~X, T?ip — > l + X, and vi°n — > l + X processes, respectively. The 
thick lines are for I = r and the thin lines are for I = \x. 



In FigEJ we plot the total cross sections of the DIS (IS) process for v T X and v T N 
scatterings by thick lines. Those of the v^N and v^N scattering processes are shown by 

^Slightly more complicated rescaling low has been examined by Albright and Jarlskog|12|. 
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thin lines for comparison. Those curves are obtained by using the parton distribution 
function(PDF)'s of Martin et al. [21]. The results are similar to the RES case, [i pro- 
duction cross sections grow rapidly from low E v , and the r production cross sections 
grow mildly from around E v = 5GeV. These results are consistent with the calculations 
of Kretzer and Reno^T], which include the NLO corrections. 

Uncertainties in the total cross section due to the modification of the structure 
function W\ (Eq. ()55)l ) and in the contribution from the Q 2 < lGeV 2 region are found 
to be rather small. A more quantitative study of the uncertainty in the theoretical 
predictions will be reported elsewhere. 

In Fig0 we show the total cross section of all the r ± production process for the 
isoscalar target. The cross sections normalized to the neutrino energy are plotted 
against the neutrino energy. The left figure is for r~ production and the right figure 
is for r + production. We find that at medium neutrino energies, the QE contribution 
dominates the total cross section near the threshold, and the sum of the QE and RES 
cross sections are significant throughout the energy range of the future neutrino oscilla- 
tion experiments. Significance of the QE and RES contribution is more pronounced for 
the v T N — > t + X reaction shown in the right-hand figure, where the DIS contribution 
starts dominating the total cross section only above E v = lOGeV. Those trends agree 
with the earlier results of Paschos and Yu|10|. 




Figure 5: The neutrino energy dependence of the total cross section of r~ (left) and r + 
(right) productions off the isoscalar target, normalized by incoming neutrino energy. 
The contributions from QE, RES and DIS processes are shown by dot-dashed, dashed 
and long-dashed lines, respectively, and their sums are shown by thick solid lines. 
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6 Polarization of the produced r 



In this section, we show the spin polarization vector of the produced r lepton as a 
function of its energy E T and the scattering angle 9 in the laboratory frame. We show 
our results for two arbitrarily fixed neutrino and antineutrino energies, E v = lOGeV 
and 20GeV, for isoscalar targets. 




E, (GeV) 



E (GeV) 



E, (GeV) 



Figure 6: Production cross section and the r polarization of the process v T N — ► t~X 
at E v = lOGeV. E T dependence of the differential cross section (top), the degree of 
polarization P (middle) and the polar component of the normalized polarization vector 
cos8p (bottom) are shown along the laboratory frame scattering angle 6 = 0° (left), 
5° (center) and 10° (right), respectively. The Histograms in the top figures and the 
circles in the middle and the bottom figures represent QE process, solid lines show RES 
process, and the dashed lines are for DIS process. The r~ spin in the r~ rest frame is 
s— ^(— sin 9p, 0, cos#p). 



Fig El summarizes our results for the v T N — > t~X process at E v = lOGeV. The 
top three figures show the double differential cross section, Eq.(j23}, as a function of 
E T , at 9 = 0° (left figures), 5° (center figures) and 10° (right figures). The DIS (IS) 
contributions are shown by dashed lines, and the RES and QE contributions are shown 
by the solid lines. The area of the histogram for the QE process is normalized to the 
cross section. A set of three middle figures give the degree of polarization, P of Eq. (|TB|) . 
as functions of E T . In the bottom three figures, we show the E T dependence of the 
polarization direction, Eq.([lb|). by using cos^p = s z /(P/2). This suffices to determine 



14 



the polarization direction because s x turns out to be always negative, s x = —^-smOp, 
and s y = (ipp = tt). It should be noted that all the 9 figures have common horizontal 
scale. The overall phase space of the E v = lOGeV experiment in the laboratory frame 
has been shown in Fig^right). 

The differential cross sections are obtained from Eq.([24p. According to the phase- 
space plot of Figfljright), along a fixed laboratory scattering angle 9, there are two 
E T 's at which QE and RES reactions can take place. The top figures of FigEJshow us 
that the cross sections in the lower E T sides are negligible. This is because of the form 
factor suppression which is significant already at E v = lOGeV. The QE and RES cross 
sections are large at forward scattering angles, and the DIS contribution become more 
significant at large scattering angles, though the cross section gets smaller. In order to 
examine the transition between the A resonance production (RES) process and the DIS 
process, we show our predictions for RES up to W < 1.6GeV and those for DIS from 
W > 1.4GeV, allowing for the overlap. Although there is no strong reason to expect 
smooth transition, we find the tendency that our predictions for the production are 
relatively smoothly changing in the transition region. 

The degree of polarization P and the polar angle 9p are defined in Eq. (jl6j) . The pro- 
duced t~ is almost fully polarized except at the very small scattering angle. As for the 
angle of the polarization vector, the high energy r~ is almost left-handed (cos^p = — 1). 
On the other hand, the spin of low energy r~ turns around. The azimuthal angle ipp 
takes tt at all energies, which means that the spin vector points to the direction of the 
initial neutrino momentum axis. 

In order to understand the above features, it is useful to consider the polarization of 
t~ in the center of mass (CM) frame of the scattering particles. Let us consider the DIS 
process in the vq CM frame, since the vq scattering is dominant in the u T N — > t~X 
process. In this frame, produced r _ is fully left-handed polarized at all scattering 
angles. This is because the initial v T and q id or s quarks) are both left-handed, and 
hence angular momentum along the initial momentum direction is zero, while in the 
final state, the produced u quark is left-handed and hence only the left-handed t~ is 
allowed by the angular momentum conservation. This selection rule is violated slightly 
when a charm quark is produced in the final state and because of gluon radiation at 
higher orders of QCD perturbation theory. The r~ polarization in the laboratory frame 
is then obtained by the Lorentz boost. In the QE and RES processes, situations are 
almost the same as in the DIS process. In the CM frame of vN collisions, the r~ 
lepton produced by the QE or RES process is almost left-handed at all angles, for the 
CM energy of y/2ME u + M 2 « 4.4GeV for E v = lOGeV, for our parametrizations of 
the transition form factors. High energy r~'s in the laboratory frame have left-handed 
polarization because those r~'s have forward scattering angles also in the CM frame. 
However, lower energy r~'s in the laboratory frame tends to have right-handed polar- 
ization because they are produced at backward angles in the CM frame. At the zero 
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scattering angle 6 = 0° of the laboratory frame, the change in the r~ momentum direc- 
tion occurs suddenly, and hence the transition from the left-handed r~ at high energies 
to the right-handed r~ at low energies is discontinuous. Since the degree of polarization 
P vanishes at this point, the polarization vector, or the density matrices are continuous. 

Figd shows the v T N — > t + X case at E v = lOGeV. The predictions of the QE 
and RES processes are quite similar to those of the v T N — > t~X process in FigJBJ 
except that the r + polarization is almost right-handed. In the DIS process, however, 
the polarization vector of r + is predicted to be quite different from the r~ case in a 
non-trivial manner. This is because the vq scattering contribution is not so small as 
compared to the vq scattering contribution. In case of the r + production process, the 
azimuthal angle of the r + polarization vector takes (fp = at all energies, which gives 
s x — -f sin 6p and s y = 0. Therefore the r + spin vector points away from the initial 
neutrino beam axis, contrary to the r~ spin case. 




23456789 10 11 23456789 10 23456789 10 

E, (GeV) E t (GeV) E T (GeV) 

Figure 7: The same as Fig® but for the process v T N — > r + X at E v = lOGeV. The r" 
spin in the r + rest frame is s = ^(sin^p, 0, cos^p). 



In order to understand the difference between the r + and t~ spin polarization pre- 
dictions in FiglHland FiglZl we show in Fig|S]the E T dependence of the differential cross 
section at 9 = 0° for the u T N — > t~X process (a) and for the v T N — > t + X process 
(b). The contributions from the v T q or v T q scattering process are shown by dashed 
lines, those from the v T q or v T q scattering process are shown by dash-dotted lines, and 
their sum by solid lines. It is clear that the vq scattering contribution dominates the 
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v T N — > r X process, whereas for the z/j-iV — > t + X process, both vq or vq scattering 
contribution are significant. 




2468 10 2468 10 
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Figure 8: The differential cross section of the DIS process at E v = lOGeV and at 9 = 0° 
in the laboratory frame for the processes v T N — > t~X (a) and v T N — > t + X (b), where 
iV is an isoscalar nuclei. The dashed and dot-dashed lines represent, respectively, the 
contributions of the neutrino-quark and neutrino-antiquark scattering subprocesses. 
The solid lines show their sum. 

Let us consider the vq and vq scattering in the parton collision CM frame. As 
for the vq scattering, the amplitude of right-handed r + production is proportional to 
(1 + cos#) and that of left-handed r + production is proportional to sin (9, where 9 is 
the scattering angle in the CM frame. On the other hand, for the vq scattering the 
produced r + has fully right-handed polarization, and the angular distribution is flat 
in the CM frame. Next, let us consider the 6 = 0° case in the laboratory frame, which 
correspond to 9 = 0° or 180° in the CM frame. Because of the (1 + cos#) 2 and sin^ 2 
distributions in the CM frame, the r + from vq scattering has fully right-handed po- 
larization and is produced only in the forward direction (9 = 0°) in the CM frame. 
Hence, all r + 's from vq scattering are right-handed along 9 = 0° in the laboratory 
frame. On the other hand, the r + 's from the vq scattering are purely right-handed 
both at 9 = 0° and 180° in the CM frame, and hence in the laboratory frame, the high 
energy r +, s have right-handed polarization, and the low energy r +, s have left-handed 
polarization. By comparing the cross section of vq and vq scattering in Fig|Hfb), we 
find that high energy r + 's are mostly produced by vq scattering, and hence they are 
almost right-handed. But as the energy decreases the contributions from vq scattering 
increases and the degree of polarization becomes lower by cancellation. In the QE and 
RES process of r + production, the mechanism is the same as the r~ production case 
but for the sign of the polarization. In the CM frame of vN scattering, r + has almost 
right-handed polarization for both QE and RES processes at all angles. Therefore in 
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the laboratory frame, high energy r + 's are right-handed, while low energy r + 's are 
left-handed because of the helicity flip by boost. 
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Figure 9: The same as Fig® but for the process v T N —* t X at E v = 20GeV. 



Let us now show our predictions at higher neutrino energies. In FigEl and Fig lTUl 
we show our predictions for u T N — > t~X and v T N — > r + X processes, respectively, at 
E v = 20GeV. The energy dependence of the differential cross sections shows clearly 
the dominance of the DIS contribution at higher energies except at 9 = 0°. FigfTUl 
shows that the QE and RES contributions are more significant in the i/N scattering, 
as compared to the vN scattering case shown in FigJUJ The relative importance of the 
QE or RES contribution to the vN scattering persists at high energies as can be seen 
from FigJ5J The degree of the r _ polarization remains high in FigEl except for the 
special angle of 9 = 0°, and its polarization direction are essentially understood by the 
boost effect, as for the E v = lOGeV case. In FiglTUl the degree of the r + polarization 
decreases at lower r + energy in the laboratory frame for the u T N — > t + X process. 
This is understood as a result of the cancellation between the vq and vq scattering 
contributions as in the E v = 10 GeV case. 

It is notable that the produced r~ has almost 100% polarization (P w 1) at all 
energies except at around 9 = 0° while its polarization direction deviates from the pure 
left-handed direction (cos#p = —1) even at relatively high r~ energies in the laboratory 
frame. On the other hand, the r + polarization deviates from 100% at relatively high 
E T while its direction is along the right-handed direction (cos#p = +1) down to half 
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the maximum energy. Those qualitative difference between the r~ polarization in the 
v r N — > r~X process and the r + polarization in the v T N — > r + X process is understood 
as a consequence of the significance of the antiquark contribution to the DIS process 
in vN scattering. 
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Figure 10: The same as Figd but for the process v T N — > t + X at E v = 20GeV. 



7 Discussion and Conclusion 

The information on the polarization of r ± produced through the v T N and v T N 
scattering is essential to identify the r production signal since the decay particle distri- 
butions depend crucially on the r spin. It is needed in long baseline neutrino oscillation 
experiments which should verify the large — > v T oscillation, and is also needed for 
the background estimation of — > v e appearance experiments which should measure 
the small mixing angle of Vg-v^ oscillation. 

In this paper we studied in detail the spin polarization of r ± produced in v T and 
v T nucleon scattering via charged currents. Quasi-elastic scattering (QE), A reso- 
nance production (RES) and deep inelastic scattering (DIS) processes have been stud- 
ied. The three subprocesses are distinguished by the hadronic invariant mass W. 
W = M(= m N ) gives QE, M+m n < W < W cut gives RES and W > W cut gives DIS. In 
this article, we set the kinematical boundary of RES and DIS process at Wat — 1.4GeV. 

The spin density matrix of r^ 1 production has been defined and the r spin polar- 
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ization vector has been defined and parametrized in the r ± rest frame whose polar-axis 
is taken along the momentum direction of in the laboratory frame. The spin density 
matrix has been calculated for each subprocess by using the form factors for the QE 
and RES processes, and by using the parton distribution functions of Ref.j2I] for the 
DIS process. We have shown the spin polarizations of r ± as function of the energy 
and the scattering angle in the laboratory frame for u T N — > t~X and v T N — > t + X 
processes at E v = lOGeV and 20GeV. We find that the produced have high degree 
of polarization, but their spin directions deviate significantly from the massless limit 
predictions at low and moderate r energies. Qualitative feature of the predictions have 
been understood by considering the helicity amplitudes in the CM frame of the scat- 
tering particles and the effects of Lorentz boost from the CM frame to the laboratory 
frame. 

Finally, we summarize our findings in Fig^J and Fig|T2l In FiglTTl we show the 
polarization vector s of r~ for the v T N — > r~X process at E u = lOGeV on the p T cosO- 
p T sin 6 plane, where p T and 9 are the produced r momentum and the scattering angle 
in the laboratory frame. The length of each arrow gives the degree of polarization 
(0 < P < 1) at each phase-space point and its orientation gives the spin direction in 
the r~ rest frame. The differential cross section is described as a contour map, where 
only the DIS cross section is plotted to avoid too much complexity. The outer line 
gives the kinematical boundary, along which the QE process occurs. Fig^2 is a more 
visual version of the information given in Fig El Fig ^] gives the r + polarization for the 
v T N — > r + X process at E v = lOGeV, compiling the cross section and the polarization 
information of Fig|7| 

Before closing our discussion, we point out some uncertainties in our calculation. 
One is the uncertainty at small Q 2 region (Q 2 < lGeV 2 ) in our DIS calculation. In 
this paper, we used an extrapolation of the parton model calculation in this region by 
freezing the PDF's below their validity region. Because the parton model must break 
down in this region, and because our estimation of the cross section in this region is not 
small, a more careful treatment, e.g. by using the structure function data is needed. 
Another is the uncertainty in the pseudo-scalar form factors, F p (q 2 ) for the QE and 
Cg (g 2 ) for the RES processes, which are not known enough so far. Because of the large 
r mass and because of the spin-flip nature of those form factors, they can affect the 
predictions of the r ± polarization significantly. QCD higher-order corrections should 
affect the polarization in the DIS region. We plan to study those uncertainties 
elsewhere. 

We hope that this work will be useful in detecting the r appearance signal in long 
baseline neutrino oscillation experiments, and that it will also be useful in understand- 
ing the — > l^vv background for the — > v e appearance experiments. 
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Figure 11: The contour map of the DIS cross section in the plane of p T cos6 and 
p T sin6 for the v T N — > r~X process at 22^ = lOGeV in the laboratory frame. The 
kinematical boundary is shown by the thick grey curve. The QE process contributes 
along the boundary, and the RES process contributes just inside of the boundary. The 
polarization are shown by the arrows. The length of the arrows give the degree of 
polarization, and the direction of arrows give that of the r~ spin in the r~ rest frame. 
The size of the 100% polarization (P = 1) arrow is shown as a reference. The arrows 
are shown along the laboratory scattering angles, 9 = 0°, 2.5°, 5°, 7.5°, and 10°, as 
well as along the kinematical boundary. 
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Figure 12: The same as FigHT} but for u T N — > t + X case. 
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